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INTRODUCTION 


This report is intended to summarize progress made on NASA Grant NAG-1- 
392 during the period August 19, 1984 to February 18, 1985. The overall goal 
of the project is to gain a better understanding of the transient behiavor of 
heat pipes operating under both normal and adverse conditions. Normal 
operation refers to cases where the capillary structure remains fully 
wetted. Adverse operation occurs when drying, re-wetting, choking, non- 
continuum flow, freezing, thawing etc. occur within the heat pipe. 

The grant was extended for one year on January 1, 1985 and the work was 
redirected towards developing the capability to predict operational behavior 
of liquid metal heat pipes used for cooling aerodynamic structures. Of 
particular interest is the startup of such heat pipes from an initially frozen 
state such as might occur during re-entry of a spac.^ vehicle into the earth's 
atmosphere or during flight of hyptrsonic aircraft. 

PROGRAM WORK STATEMENT 

The program work statement has been changed several times during the past 
two years. The original long range statement included the following tasks. 

I. Compile governing transient equations for heat pipe operation. 

II. Model adverse conditions and select non-dimensional groups which 
can be used to define regimes where adverse operation occurs. 

III. Incorporate models for adverse operations into solution 
techniques. This will involve using thermal property subroutines 
to predict properties for use in conduction equations. 

IV. Assess adequacy of existing finite element computational schemes to 
handle the models developed. 
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V. 


Develop finite element representation of heat pipe for normal 
operation. 

VI. Include non-dimensional groups which define adverse operating 

regimes in computational scheme. 

VII. Select test case to demonstrate techniques. 

VIII. Compute operating parameters for test case. 

IX. Test computed values against experimental results for test case. 

X. Perform parametric studies for system to the extent possible. 

XI. Develop, to the extent possible with given time and resources, 

simplified correlation equations or design procedures. 

About one year ago the following tasks were added. 

I. Compile a Space Station thermal control candidate technology data 

base. 

II. Develop methodology and algorithms to evaluate candidate 

technologies. 

III. Demonstrate methodology and algorithms through computer-aided 

assessment example using data base. 

The grant was funded for an additional year starting January 1, 1985 and 
the emphasis was changed to the following tasks. 

I. To investigate the physics of melting and freezing of the working 
fluid in liquid metal heat pipes, and of other performarce limitations which 
are presently not well understood, and to develop equations, adequate for 
design, to predict associated performance limits. 

II. To evaluate the Langley finite difference heat pipe design computer 
program and the equations upon which it is based and suggest improvements or 
additions to upgrade the program. 
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III. To work with the aerospace contractor who wins the heat pipe study 
contract to upgrade the contractor's heat pipe analysis capability. 

IV. To work with NASA to define the required instrumentation for the 
heat pipe model which will be built under the heat pipe study contract and to 
develop a test plan for the heat pipe. 

While these additions and changes in emphasis in the work statement may 
at first appear to be excessive, the overall goal has not changed. That goal 
is to develop predictive techniques for transient heat pipe operation. 

EVALUATION OF SPACE STATION THERMAL CONTROL CANDIDATES 

A data base is being compiled for candidate technologies for thermal 
acquisition, transfer and rejection for the Space Station. Representative 
data were presented in the last progress report. Additional data have been 

compiled, consolidated and assessed, and typical results are presented in 
later sections of this report. 

During the period covered by this report a computer-aided assessment 
program was obtained from John B. Hall, Jr. of NASA/Langley Research Center. 
This program, written in FORTRAN for a Prime Computer, has been installed on 
the CDC Cyber 180 at Georgia Tech. The installation required a significant 
number program modifications and was completed in January 1985. Details of 
the modifications can be provided upon request. The algorithms useo for the 
assessment of candidate technologies have been identified and are currently 
being evaluated. Examples of the evaluations are discussed in the following 
paragraphs. 

An assessment of candidate data for two thermal acquisition candidates, 
the conductive cold plate and the two-phase cold plate, has been initiated. 
The present evaluation is based upon the data in Ref. [1] wherein the design 
of the internal thermal bus of a space platform is considered. Included in 


TABLE 1 


Thermal Control Subsystem Heat Loads (after Ref. [1] ) 


Module 


Heat Load, kW 


Metabolic 

Loop 

Equi pment 
Loop 

Subsystem 

Total 

Logistics 

2.36 

1.0 

3.36 

Habitat 

2.36 

5.0 

7.36 

MBA 

2.36 

6.0 

8.36 

Lab. No. 1 

2.36 

10.0 

12.36 

Lab. No. 2 

2.36 

30.0 

32.36 


TABLE 2 

Thermal Control Subsystem Weights and Pump Power Requirements 
for Conductive Cold Plate (after Ref. [1] ) 


TCS Weight, lb 


Pump Power, W 


Module 

Metabolic 

Loop 1 

Equipment 

Loop 2 

Subsystem 

Total 

Metabol 1 c 
Loop 

Equi pment 
Loop 

Subsystem 

Total 

Logistics 

69.4 

"84.0 




45.6 

Habitat 

69.4 

226.3 

295.7 

21.8 

102.0 

123.8 

MBA 

69.4 

276.0 

345.4 

21.8 

123.0 

144.8 

Lab. No. 1 

69.4 

452.0 

521.4 

21.8 

194.0 

215.8 

Lab. No. 2 

69.4 

1213.5 

1282.9 

21.8 

566.0 

587.8 


1 . 


2 . 


Metabolic loops employ air/water heat exchangers to cool cabin air. 
Equipment loops employ conductive cold plates for equipment cooling. 





the platform are five pressurized modules whose heat loads are summarized In 
Table 1. The metabolic loops for each module are used to provide temperature 
and humidity control for the crew metabolic loads by means of air/water heat 
exchangers. Each module also has cold plates to provide equipment cooling. 

CONOUCTIVE COLD PLATE SYSTEM 

The following assumptions were made for the conductive cold plate systems 
analyzed in Ref. [1]: 

• The average heat flux density for the conductive cold plate is 
270 W/ft 2 . 

t The cold plate density is 5.3 lb/ft 2 . 

• The pumping power penalty for the cold plate Is 350 lb/kW. 

• Water Is the working fluid. 

• Temperatures at the cold plates are limited to 20 5°C. 

The total pumping power was determined from a design analysis considering 
pressure drops throughout the system. Table 2 summarizes the results of the 
analysis with respect to the weight and pumping power for each of the five 
pressurized modules. These data and the corresponding loads from Table 1 are 
presented graphically in Figures 1 and 2. Included in these figures are least 
squares curve fits which can be incorporated into the computer-aided 
assessment program. Other candidate data derived from Ref. [1] include the 
fol 1 owl ng: 

• Cold plate volume ( f t ^ ) = 0.11 * Heat Load (kW) 

t Cold plate heat transfer surface area = 3.7 ft 2 /kW heat load 

TWO-PHASE COLD PLATE SYSTEM 

The following assumptions were made for the two-phase cold plate systems 
analyzed in Ref. [1]: 

• The average heat flux density for the two-phase cold plate is 
600 W/ft 2 . 
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Subsystem Weight, W (lb) 







9 The cold plate density Is 11.5 lb/ft 2 , 

e The pumping power penalty for the cold plate Is 350 Ib/kW. 

• Water Is the working fluid and complete vaporization occurs In 
the cold plate. 

e Temperatures at the cold plates are limited to 20 5°C. 

Table 3 summarizes the results of the design analysis from Ref. [11 with 
respect to the weight and pumping power for each of the five modules. These 
data and the corresponding loads are presented graphically In Figures 3 and 
4. Included In these figures are least squares curve fits which can be 
Incorporated Into the computer-aided assessment program. Other data derived 
from Ref. [1] Include the following: 

• Cold plate volume (ft^) * 0.175 * Heat Load (kW) 

• Cold plate heat transfer surface area ■ 1.67 ft 2 /kW heat load 

COMPUTER IMPLEMENTATION OF CANO 10 ATE EVALUATION ALGORITHM 

The candidate data bases are currently set up for a candidate power 
rating of 50 kW so that all data are based upon this rating regardless of 
whether the system Is used for thermal acquisition, transport, or rejection. 
Mission data are used to define the actual power rating, and candidate data 
are scaled up or down depending upon whether the mission power rating Is 
greater than or less than 50 kW, respectively. 

This scneme could be retained for use with new candidate evaluation 
algorithms. Defining a load fraction, 

LFRAC ■ SIZE/50 

where SIZE Is the mission power rating, the pump power required, for example, 
could be written as 

P s ■ P s (LFRAC) 


- 8 - 


TABLE 3 


Thermal Control Subsystem Weights and Pump Power Requirements 
for Two-Phase Cold Plate (after Ref. [1]) 




Ter Weight, lb 


Pump Power, W 



Metal bol 1c 
Loop 1 

Equipment 

Loop 2 

Subsystem 

Metal bol 1c Equipment 

Subsyste 

Module 

Total 

Loop 

Loop 

Total 


|]R1 

97.1 

■tuu 

BtRiH 

BiWQB 


Habitat 

91.2 

198.9 

290.1 

0.04 

0.17 

0.21 

MBA 

91.2 

222.2 

313.4 

0.04 

0.26 

0.30 

Lab. No. 1 

91.2 

313.9 

405.1 

0.04 

0.34 

0.38 

Lab. No. 2 

91.2 

757.5 

848.7 

0.04 

0.86 

0.90 


1. Metabolic loops employ air/water heat exchangers to cool cabin air. 

2. Equipment loops employ two-phase cold plates for equipment cooling. 
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Subsystem Weight, W (lb) 



Figure 3. Variation of Subsystem Weight with 
Heat Load - Two-Phase Cold Plate 
System . 
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1.0 



Figure 4. Variation of Subsystem Pump Power 
with Heat Load - Two-Phase Cold 
Plant System. 
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1 




where the function on the left side of this expression is a least squares 
curve fit of data from the design analysis. The pump power required for the 
actual mission size could be calculated in the computer-aided assessment 
program in a number of ways, for example: 

Scheme 1 - Add a subroutine or function subprogram which defines P s = 

P $ (LFRAC) for each candidate. This would probably be the 
least flexible means of implementation because the program 
user could not easily create the subprograms for the various 
candidates. 

Scheme 2 - Store the coefficients of a least squares polynomial curve fit 

(order <_ 2 should be sufficient) in the data base. This would 
at least require (a) user input of up to three coefficients 
used to define P s ( LFRAC ) when creation of a data base is 
desired, (b) larger array sizes for the data arrays used in 
the program (currently dimensioned to 82), (c) revision of the 
indexing scheme and dimension for the SYSRAY array in 
subroutine COMPUTE, (d) revision of the equation used to 
compute the power required in subroutine COMPUTE; for example, 

SYSRAY( INDEX+3) = DAT ARAY ( N1 ,1) + (DATARAY(N2,I) + DATARAY(N3,I )*LFRAC)*LFRAC 
where N1 , N2, N3 are the indices of the array elements 
containing the coefficients for the least squares 
polynomial. All other re iationships derived from design 
studies or analyses (i.e. relationships for weight vs. load, 
volume vs. load, heat transfer surface area vs. load, etc.) 
could be implemented in a similar manner without much 
di fficulty . 
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Liquid Film 




The computer-aided assessment program installed at Georgia Tech is 
currently being modified along the lines described in Scheme 2 above. 
Complete details of the required modifications as well as sample assessment 
results will be reported as soon as they have been completed. 

MODELING OF HEAT PIPE STARTUP FORM THE FROZEN STATE 

Liquid metal heat pipes experience operating difficulties when a portion 
or all of the working fluid is solid. Of particular interest, in the 
utilization of these devices, is startup from an initial state where the 
working fluid is totally frozen. This problem has been considered by Cotter 
[2], Ivanovskii et al. [3], Deverall et al. [4], Neal [5], Camarda [6], 
Silverstein [7] and others. However, no comprehensive experimental or 
theoretical investigation of the problem has been performed. The goals of the 
present work, therefore, are to examine the physics of the processes that 
occur during startup from a frozen state, to mathematically model the 
processes, to develop numerical techniques for evaluating the mathematical 
models, and to check the mathematical models and solution techniques against 
experimental data which will be provided by NASA. 

Currently, it is believed that the following sequence of events occurs 
during startup from a frozen state. Initially the working fluid is solid and 
the vapor density is extremely low so that molecular flow conditions prevail 
throughout the heat pipe vapor section. When thermal energy is added in the 
evaporator section melting of the working fluid begins and energy transport 
from the heated zone to the adjacent unheated zone proceeds quite slowly via 
axial conduction through the heat pipe wall and the combination of working 
fluid and wick structure. However, energy transport by the vapor is 
negl i gi bl e. 
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Since energy is continuously added at the evaporator section, after a 
time the vapor density in the heated zone is increased so that the mean free 
path of molecules becomes small compared to the diameter of the vapor 
passage. Thus, in the hot zone a continuum flow regime is established, and in 
the cold zone the vapor space is still in free molecular flow. The continuum 
flow regime grows with time until it reaches the end of heat pipe. 

After the frozen working fluid in the evaporator section is melted, 
evaporation can take place at the liquid-vapor interface and vapor can flow 
toward the condenser section. Initially vapor generated in the evaporator 
freezes uniformly on the inner surface of the frozen working fluid in the cold 
zone. For this stage, energy is transferred as latent heat owing to 
vaporization in the heated zone and condensation and freezing in the cooled 
zone. 

As long as the entire condensate layer is below the freezing temperature, 
it remains in a solid state. However, the vapor-solid interface temperature 
increases as the vapor deposition proceeds until it reaches the melting 
point. Since the central vapor temperature is above the freezing temperature, 
further condensation of vapor causes simultaneous melting and condensing at 
the interface, and a liquid layer appears near the leading edge of the frozen 
working fluid. 

During this part of startup the vapor pressure gradient is large and as 
shown in Figure 5 the vapor flow passage at the condensing zone has the shape 
of a convergi ng-di vergi ng nozzle. The velocity of the vapor is so high that 
the vapor shear acting on the interface between the liquid film and vapor is a 
major driving force of the liquid film. This liquid film is continuously 
dragged downstream from the leading edge until the liquid meets a colder 
temperature to solidify. This annular type flow is assumed to be the dominant 
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flow pattern of the liquid film existing over the condensing zone. In this 
case, failure of a heat pipe may occur because freezing of the condensed vapor 
prevents liquid return to the evaporator. 

This process continues until the frozen working fluid is completely 

melted. When the entire working fluid is finally melted, sufficient liquid is 
returned to the evaporator and normal transient operation begins. Eventually 
the heat pipe reaches near isothermal and steady state operation. 

On the basis of experimental observations, startup behavior of a frozen 
heat pipe may be divided into three distinct phases for convenience of 
analysis. 

Phase I: The vapor flow in the heat pipe is in a free molecular 

regime throughout the vapor space. 

Phase II: In the vapor space, a continuum flow regime is established 

in the heated zone and a continuum flow front moves toward 
the other end of the heat pipe. 

Phase III: Continuum flow exists over the entire heat pipe length in 

the vapor region. 

Phase I Modeling 

Heat transfer by the vapor is negligible and heat tranasfer from the 
heated zone to the unheated zone proceeds quite slowly by conduction only. 
Hence, conduction can be accepted as the sole transport mechanism. The heat 
pipe is subcooled and initially has uniform temperature, T . Continuum flow 

ao 

is considered to exist in the vapor space when the temperature of the vapor 
space is greater than that calculated by the following equation. 
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when u, p , R and D are viscosity, density, gas constant and diameter of vapor 
space, respectively. 


Heat transfer is governed by the conduction equation written in 
cylindrical coordinates as follow: 

r 5 T' + !l< K iir> ■ <•»'), S 1 . 1 ■!>.».* (2) 

where the subscripts p, i , and s denote the heat pipe shell, liquid region, 
and solid region, respectively. 

At each interface, the temperature and the heat flux are continuous. 
These boundary conditions are 


1. At r = r w (Interface between the shell and wick) 


3T ) , 

K_ = K — ^ and T D = T 
P 3r i 3r PH 


: liquid region 


3T P 9T s 

K P 3 ^ * K s TT and T P * T s 


: solid region 


T P ■ T s ■ T , 


: subcooled solid region 


where T is the initial temperature. 


2. At r = r (Liquid-solid interface) 


T = T = T (melting temperature' 
l s m J 


From an energy balance at the interface, 
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3T 36 3T 

K i F s p s h sn It + K s Jr~ 

where 5 $ is the position of the liquid-solid interface. 

3. At r = r (t) (forward edge of the solid phase temperature wave) 

oo 



This boundary condition is valid as long as r > r. . When r < r. , heat 

00 ] 00 

is transferred to the free molecular region. 


4. At z = 0 and L for all r (i.e. both ends of the heat pipe) 



5. No heat transfer occurs external to the shell at the adiabatic 
section 



6. The external surface of the shell at the evaporator and condenser 
sections can be exposed to a heat flux, convection or radiation. 
Thus 

k p TT - n(T - ' V + ° F(T - 4 ' t p 4) + q " 
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or by conduction alone 



where subscript o stands for surrounding material. However, during 
Phase I, no heat transfer occurs at the condenser section. Thus 

3T P 

- — = 0 at condenser section 

3r 

Phase II Modeling 

A. Liquid flow in the wick 

After the frozen working substance has melted in the evaporator the 
liquid-vapor meniscus has a more concave surface due to evaporation of working 
fluid. Thus a pressure gradient is set up and the liquid generated owing to 
phase change at the liquid-solid interface in the wick flows back to the 
evaporator through the wick. When the wick is no longer able to supply 

sufficient liquid, the liquid-vapor interface retreats further into the wick. 

The analytical model includes the effects of the rate of change of 
momentum, surface tension forces, frictional forces in the body of the wick as 
well as axial variation of static vapor pressure. The following assumptions 
are made concerning the liquid flow model. 

1. The wicking material is isotropic, is of constant thickness, and is 
initially saturated with a wetting liquid. 

2. The liquid flow is axi symmetric , incompressible, laminar, and 

unsteady. The liquid velocity V 1 is uniform across the wick, and has 
a z component only. 

3. The liquid-vapor interface meniscus can be characteri zed by one 
radius of curvature at axial location. 



4. The temperature of the liquid Is a function of 2 and r, and the 
kinetic energy and the effect of gravity are neglected. 

5. The vapor tha- evaporates or condenses on the liquid-vapor Interfaces 
has a velocity In the radial direction only. 

Application of the principle of conservation of mass for the fluid 
element shown in Figure 7 yields. 


3_ 

3t 


(m ) * — (rn )dz - m 
' j r 3z ' v v 


(3) 


where 

% * p AV "’ll * p A dZ * m v = u o p » PdZ > A c * e " (r w - r i>- p * 2 ” r i 


and u Q is the velocity of the vapor. Substitution of these expressions 

• • 

for m , m , m , A , and P into equation (3) yields 

** V C 



3r 1 < r w 2 - r i> 3V , 


a 3z 


7F 


dz 


U p 
0 M V 

ep. 


(4) 


When the entire wick is saturated, equation (4) becomes 

3V t 2 ” r i % p v 



(5) 


The momentum equation in the axial direction for the fluid element in 
Figure 7 is obtained from Newton's second law with Laplace-Young's equation 
for capillary forces and Darcy's law for friction force. 


L F 


z 


= M 


z-dz 


M + M 

Z V 
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Figure 7. Differential Element of Wick Material for Mass Balance and 
Manenturr Ecuation. 



Differential Element of Wick Material for Energy Equation. 


- 2 


Figure 8. 



or 


Jr * v t IT - r . [ TT * 2 ° h < I » 


3V 


aP. 


4 32 o u 32 


2r, 


3 r\ 


-1 /. 2o * 1 . e #4 xw 

, r 5 .it (p v * r ) ir + r < r )v t 

P t (r w - r,) 1 ' 


( 6 ) 


where Is the permeability of the wick and o Is liquid-vapor surface 
tension. 

The principle of conservation of energy applied to a differential wick 
element as shown In Figure 8 yields 


It <Vi> ■ It (A c K e > dz + If (2 ’ rK e Jr dz)dr 


+ H (m t h t )dZ 


(7) 


or 


3T l 

(tB t c v> ir 


a 3 "!’. i a 3 T 3 T 

Tz^eTT ) r fr ^ rK e 7r^ ' + <« ( C P )V t aT 


+ £ iL ( r y p ) 
r 32 V 2 4 ; 


( 8 ) 


where K g is the thermal conductivity of the liquid-wick combination. 


Boundary conditio: i are 


1. at Z * 0 


3T 4 

V ■ 0 and — -=• ■ 0 
4 32 
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2. at r « r w (interface between liquid and tube shell) 


at aT. 


and T « T 


p Jr e ar P « 


3. at r « r (liquid-vapor Interface) 


Boundary conditions will be stated In section B. 


4. at z * z.j (liquid-solid interface In wick) 


r v * r 1 


p s 3Z 1 

V « — — L 

i p t at 


- K 


!L 

e 3z 


h 

P S h St 3t 


- K 


s 3Z 


(mass balance) 


(energy balance) 


T 

l 


T 

m 


(melting temperature) 


B. Vapor flow 

It is assumed that the vapor flow is ax i symmetric, compressible, and 
laminar. The vapor flow can be described by the Navi er-Stokes equations for a 
compressible fluid by neglecting body forces and assuming constant viscosity. 


Continuity equation 


_3p 

3t 


+ 7f?' rou 



« 0 


(9) 
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Momentum equation for the r-dlrection 


r SU * 3U * Y/ 1 « 3 P 4 . 4 ,3 r 1 9 

^it “ir V iz* ir ( IF C r 7r (ru)]) 


. 1 3 , 3V , S 2 U 

+ ? u Tr < U > * u ; Z T 


( 10 ) 


and the z-di rent Ion 


o[ 3V* u iV + v S -i]-.^ +u [^ + iH t l!| ] 

PL at 3 r az J 3z u L 3f 2 r 3 r az 2 J 


Energy equation 


+ 7 “iz £?!?<*’“> 


(in 


3T 


3T 


3T. 


(pC, v [ ir * u i?" + v it ] * ? If (r K w tt 


37 


3 3T v 

+ — (K — - ) 
3Z v v 3Z ' 


( 12 ) 


Whether there Is a liquid film or not, two different boundary conditions exist 
in the condensing region. When no liquid film exists 

1 . at z * 0 

37 

u * 0, V * 0, — - = 0 

3 Z 

2 , at r » 0 (at center of heat pipe) 
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3V 3T v 

u = °* If • °- ir 1 ' 0 


3. at r = r y and 0 < z < L g (liquid-vapor interface at evaporator) 


V = 0, T = T = T . 

i v sat 


1 


3 i i 

tv 


k- i" ( r > 

3 


where and are the pressure and temperature at the triple point. 


U • - U„(z,t) = - [K ~r ) - K il )„] 


»v\* e 3r 


3r 'v 


4. At liquid-vapor interface in adiabatic section 


3T t 3T v 

u = 0, v = 0, K — ^ = K — - 
e 3r v 3r 


5. at r = r.j - 5 $ and z^ < z <; z f (solid-vapor interface) 


v = 0 


p 35 . 3T 3T 

u - u o (z,t) = ( — • ) gt” ’ h ys p v ^ K s ~ " K v 3 r~ ^ 


T s - T v - *sub 


1 


1 R ln / 3 . 

TJ-K^ ,n ( r> 


6. At the end of continuum flow (z = Zf) 


3T 

u = 0, V = 0, — - = 0 
0 z 
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If a liquid layer exists on the solid layer, other boundary conditions are 
needed. These conditions will be stated in the next section, 

C. Liquid film 

When the temperature at the solid-vapor interface reaches the melting 
temperature, a thin liquid film is formed over the solid by simultaneous 
condensing of vapor and melting of solid. The following assumptions are made 
for mathematical formulation of the model shown in Figure 6. 

1. A continuous laminar film exists. 

2. Fluid property variations and gravity effects are neglected. 

3. The transient and convection terms are neglected in the momentum and 
energy equation. 

Continuity equation 

H ♦ It < f ’ V r > 

J r i - s 

Momentum equation 



3 3V f 
L- ( r — L) = o 

3r u ar ; 


(14) 


Energy equation 


a_ 

ar 



= 0 


(15) 


Boundary conditions 


1. at r * r. - 5 (liquid - solid interface) 
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6 = 6 s = 0 for al 1 time 
3. at z = Zf 


6 =0 for al 1 time 

s 

4. at r = r. - 6 (liquid-vapor interface) 

u Tr )y * u I? )f ( cont ^ nuit ^ shear stress) 


V y = Vf (continuity of velocity) 

p (V + u ) = m 
M v ' v 3z v' c 

h v£ m c = K fr )f ’ K It )v ( interface energy balance) 


0. Solid layer 

The governing equation is the same conduction equation that describes the 
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temperature distribution of the frozen substance in Phase I. Also, 
conditions have been stated in the previous section. 


boundary 
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